Nonperturbative renormalization group approach to Lifshitz critical behaviour 
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The behaviour of a d-dimensional vectorial N = 3 model at a m-axial Lifshitz critical point is 
investigated by means of a nonperturbative renormalization group approach that is free of the huge 
technical difficulties that plague the perturbative approaches and limit their computations to the 
lowest orders. In particular being systematically improvable, our approach allows us to control the 
convergence of successive approximations and thus to get reliable physical quantities in d = 3. 
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Lifshitz critical behaviour (LCB) [1] (see also[2-4]) oc- 
curs when a disordered phase encounters both a homoge- 
neous ordered phase and a spatially modulated ordered 
phase with a modulation wave- vector q mod ^ 0. In the 
general case the vector q mod spans a m-dimensional sub- 
space of the d-dimcnsional space with < m < d. For a 
A-component order parameter the universal behaviour at 
criticality is completely determined by the set (m, d, A). 
LCB has been proposed to occur in many systems in- 
cluding magnetic models (notably the ANNNI model [5]), 
liquid crystals, microcmulsions, polymer mixtures, ferro- 
electrics, high-T c superconductors, see [4] for a review. In 
the domain of magnetic materials there has been a grow- 
ing activity in the search for LCB behaviour. A clear-cut 
LCB has been found in manganese phosphide (MnP) [6] 
and, possibly, in the ternary uranium silicidc (UPD2Si2) 
[7] . One can thus expect accurate determinations of the 
critical quantities from experiments in a near future. 

From the theoretical point of view the simplest model 
displaying LCB can be obtained by generalizing the 
Hamiltonian, or action, relevant to study the usual vec- 
torial ferromagnetic-paramagnetic phase transition. Let 
us consider a A-component vector field <$>{x) in a d- 
dimensional space. The coordinates x are decomposed 
into a parallel component xu G ]R m and an orthogonal 
component x±_ G TR d ~ m , i.e. x = (xu,xx)- The action 
allowing a LCB writes: 



*x ± d r ' 




(i) 

where d\\ and d± stand for the derivatives in the cor- 
responding directions. The coupling constants Zu , Zj_ 
and u are supposed to be always positive while po an d K 
are allowed to change sign. The coupling K stands for a 
magnetization occuring in the - homogeneous - ordered 
phase. From a mean- field analysis one observes that, for 
po > 0, when the coefficient r = — un in front of <p 2 varies 
from a positive to a negative value the system undergoes 



a phase transition from a disordered to a homogeneous or- 
dered phase while for p$ < &, transition occurs, for some 
t = r c , from a disordered to a modulated ordered phase. 
The two transition lines join at the Lifshitz point which, 
within a mean-field analysis, is located at r = po = 0. 

The salient property characterizing LCB is that of 
anisotropic scale invariance (ASI). Indeed, at the Lifshitz 
point, because of the absence of (d^cp) 2 term, the scaling 
dimensions in the _L and || directions differ. In particular 
the two-point correlation functions scale as [4]: 



r( 2 >(q ± = 0, q | -+0) 



(2) 



which define the two anomalous scaling dimensions r]± 
and 7711 . On the other hand, for a generic scaling operator, 
one expects the following asymptotic behaviour under 
a scale transformation: 0(sq_i_, s z q||) ~ s~ A 0(q±, q||) 
when s — > where A is the scaling dimension associated 
to the operator O, z being the anisotropy critical expo- 
nent. In particular, for the two-point function one has: 
r( 2 )(sq ± ,s 2 q||) ~ s 2 -'^r( 2 )(q ± ,q||) when s -» 0. This 
behaviour, together with Eq.(2), provides the relation: 



2-y± 
4 - m 



(3) 



Finally two critical exponents and v±_ characterize the 
behaviour of the correlation lengths near criticality: 
£|| oc t~^ii and £j_ °c t~ u± with vu = zu±. 

ASI occurs in many contexts: equilibrium critical phe- 
nomena of anisotropic systems, like those described by 
action (1) at a Lifshitz point or, e.g. in the crumplc- 
to-tubule transition in anisotropic membranes [8, 9], as 
well as in dynamical critical phenomena at and away 
from equilibrium (see [4]). In quantum field theory an 
intensive activity has been developed towards theories 
for which Lorentz invariance is broken at high energy 
by high order derivative terms in the spatial directions 
(see [10] for a review). In these "Lifshitz-type theories" 
the presence of anisotropy between temporal and spa- 
tial directions drastically improve the UV behaviour and 
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renormalizability properties. These ideas have been fur- 
ther extended towards anisotropic scale invariant grav- 
ity [11] and cosmology [12]. Finally a theory of local 
scale invariance (LSI) has been introduced [13] both for 
equilibrium and out of equilibrium phenomena leading 
to conjecture exact expressions for the two-point corre- 
lators of anisotropic systems. While Monte Carlo results 
[14] have been claimed to agree with these predictions, 
a very recent work [15] based, for a part, on a pertur- 
bative approach at two-loop order in the vicinity of the 
upper critical dimension seems to invalidate LSI and its 
relevance to describe the physics of anisotropic scale in- 
variant systems. 

In this context it is clear that an efficient and system- 
atically improvable approach of anisotropic systems, and 
in particular of LCB, is needed. From this point of view 
one has to emphasize that the available, perturbative, 
techniques arc especially in trouble. Let us start with 
the weak-coupling e-expansion. A first problem is that 
going from isotropic to anisotropic systems shifts the up- 
per critical dimension from d uc = 4 to d uc = 4 + m/2. 
This means that even for the minimal non trivial value of 
m, equal to one, the e-expansion implies to deal with the 
large value e = 3/2 when computing the critical prop- 
erties in d — 3. Assuming that the series obtained are 
Borel-summable, what is not guaranteed, getting reliable 
physical quantities thus implies to compute, at least, up 
to four or five-loop order. But then one faces a second 
and important problem. As emphasized in [4, 16, 17] 
the real-space free propagator used to perform the e- 
cxpansion takes a very complicated form, known as Fox- 
Wright generalized hypergeometric functions, that leads 
to enormous [18] technical difficulties. This explains why 
the weak-coupling e-cxpansion results have been very 
controversial during a long time [17, 19] and that it took 
almost twenty years to fill the gap between early one-loop 
order results [1] and the complete two-loop order compu- 
tation [4, 16, 17, 20]. For this reason it is extremely 
unlikely that higher order contributions will be obtained 
in a near future. Similar difficulties occur within a large- 
TV approach (see [21] for a review) and it is only very 
recently [18] that consistency between this large- N ap- 
proach at order 0(1/ N) and the weak-coupling expan- 
sion at two-loop order has been firmly established. Fi- 
nally note that it is also possible, in principle, to investi- 
gate LCB by means of a low-temperature approach in the 
vicinity of the lower critical dimension which, for N > 1 
components system, is given by di c = 2 + m/2. This has 
been done at one-loop order by Sak and Grest [22] . How- 
ever, as in the O(N) model, the series obtained within a 
low-temperature approach are generally suspected to be 
non-Borel summable and thus of no practical use. 

We investigate here the LCB by means of a nonper- 
turbative renormalization group (NPRG) approach. Our 
computation is based on the concept of running effective 
action [23] (see [24-27] for reviews), [</>], a functional of 



the TV-component vector field </>(x) that describes the ef- 
fective physics at a coarse grained scale k. Technically the 
index k stands for a running scale that separates the high- 
momentum modes, with q > k, from the low-momentum 
ones, with q < k and Tk[(f>] represents a coarse grained 
free energy where only fluctuations with momenta q > k 
have been integrated out. The running of k towards k = 
thus corresponds to gradually integrate over all fluctua- 
tions. The fc-dependence, RG flow, of Tk is provided by 
an exact - albeit one- loop - evolution equation [23]: 

where t = ln/c/A, A being some microscopic, lattice, 
scale. The trace in (4) involves a d-dimcnsional mo- 
mentum integral over a momentum q as well as a sum- 
mation over vectorial indices. The function Rk(q) real- 
izes the split between low- and high-momentum degrees 
of freedom while 1^ represents the second functional 
derivative of Tk with respect to <f>, i.e. the inverse field- 
dependent propagator. Considering Tk in its full gener- 
ality Eq.(4) provides an exact RG flow for the coupling 
constants associated to any power of <fi and of its deriva- 
tives. 

There are several major advantages in using Eq.(4). 
First, one deals with a one-loop equation while the com- 
putations are naturally performed in momentum space. 
In this way we avoid all the technical difficulties en- 
countered within the perturbative approaches that are 
associated to the multi-loop structure and the complex- 
ity of the propagator in real space. Second, the equation 
being nonperturbative in the coupling constants enter- 
ing in the effective action (the 4 -like coupling constant 
it, the temperature T ~ 1/k as the parameter 1/N) the 
approach overcomes a major problem of the perturba- 
tive theory: the need to resum, if possible, the pertur- 
bative rcnormalizcd series. Third this technique is sys- 
tematically improvable without conceptual or technical 
difficulty. Let us develop these last two points. Eq.(4), 
although exact is not exactly solvable. One thus has to 
consider truncations of r *;[</>] and, thus, approximations 
of Eq.(4). Different kinds of approximations are allowed 
that keep untouched the nonperturbative and one-loop 
character of the equation. A very useful and efficient ap- 
proximation is based on an expansion of Tk [<p] in powers 
of both fields and field derivatives. A derivative expan- 
sion is particularly justified to investigate critical phe- 
nomena whose physics is dominated by low momenta and 
thus by low powers of the field derivatives. A field ex- 
pansion is justified by both its general character and its 
rapid convergence [28] apparently without need of resum- 
mation procedure, as attested by several studies involv- 
ing Ising model [29], frustrated magnets [26], randomly 
dilute Ising model [30] and other [25]. Moreover in the 
context of anisotropic systems a field expansion is par- 
ticularly suitable since the physical dimension 3 is close 
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to the lower critical dimension of (vectorial) anisotropic 
systems di c = 2 + m/2 in the vicinity of which, the RG 
flow (4), together with the ansatz (1), is one-loop ex- 
act. One can thus expect the ansatz (f ), or a bit more 
sophisticated ansatz, to provide very sensible results in 
d = 3. In this article we provide the RG equations for 
the coupling constants entering in the action (1) while 
we have computed with powers of the field up to order 
cf> 12 . The validity of this approach is then checked by 
studying both the cut-off independence of the physical 
quantities and their behaviour when the field content is 
enriched. We show in particular that converged critical 
quantities are obtained using a limited number of powers 
of the field. Note that a NPRG approach of LCB in the 
Ising case has already been performed in [31] using a full 
but local potential approach of the Polchinski equation, 
thus neglecting the anomalous dimension. We investigate 
here the behaviour of the vectorial N = 3 case providing 
both the critical exponents v\\ and v± together with the 
anomalous scaling dimensions r]\\ and rj±. 

The flow equations for the coupling constants 
K, u, po entering in (1) arc obtained, as usual 
[25, 26] by appropriate functional derivatives of the 
RG equation Eq.(4). One defines the dimension- 
less quantities: 7t = Z 1 ^ ™)/ 4 ^m/4^(m+4 2rf )/ 2 ^ 

u = Z ( r- 8)/4 Z- m/4 k { l d - m - 8]/2 u and p = 

^ i<2 1/2 1 

Z ±_ kj_ po and their flow writes, with t =ln 

fcx/A: 

d{K = — (d — m + z(m + 7711 — 4)) k 
+(N - l)f ± 2 + 37° 2 

diu = (d — m + z(m + 27711 — $)) u 
+2u 2 f(N - l)t ±A + 9lJ i4 ) 

dtPo = z(t?|| - 2)p 

(5) 

while the running anomalous dimensions r\^ and r)\\ are 
given by: 

1 f-fi 70 \ 1 f-fi -fi \ 

V±= =(^,2 + *|| 1 2j-^(^,0-'||,0j 

?? " = 6zm(l + 2) ¥¥ ( (m + 2) ^ ~ ^) 

-2(m + 2)u 3 n 3 (s]_ fi + Sl 6 ) - \ (t!, 2 - T\, 2 ) 

+9u k (t\ a + rj i4 ) - 8u 2 k 2 (t% - tJ i6 ) 

+4U 3 k 3 (t% + tJ )8 ) + 6u 2 k 2 (Fi i2 - uy 

-12U 3 k 3 (ul A + ul A ^ 

_ a _ a _ a _ a _ a (6) 
where, in Eqs.(5) and (6), l a bl M™ <b , S^ b , T" & , U^ b are 



dimensionless "threshold functions" (see [25]) l" b , M" b , 
Sl b , T« b , U- b that are given by : 

"a r K b q 2a F(q,|) 

9tJ [P(q||)+m2q|] 

where K b = -r[ 7b ](47r)( m - d )/ 2 /2, 7p = (m + p - d)/2, 
P (<1\\) = z \\ q| + Poqf + Rk±(<l\\), m 2 = 4uk, m_L = 
and where d/dt only acts on R^ ± . In Eq.(7) the function 
F(q,|) is given by 1, (dP/dqjj) 2 , (dP/dqj) 3 , {dP/dq\) A 

and (d 2 P/d(q|) 2 ) 2 for /, M, S, T and U respectively. 
These threshold functions encode the nonperturbativc 
content of the approach since, as it is clear from (7), 
they are nonpolynomial functions of the coupling con- 
stants u and k entering in the squared "mass" m|. Note 
that the threshold functions only depend on qii since the 
momenta q^ have been integrated out. 

Let us discuss the RG equations Eqs.(5) and (6). The 
one- loop structure of Eq.(4) together with action (1) al- 
low, when the Eqs.(5) and (6) are expanded in powers 
of the suitable coupling constant, to recover all the one- 
loop results obtained perturbatively. The weak-coupling 
results obtained in the vicinity of the upper critical di- 
mension [1] are easily recovered by performing an expan- 
sion in powers of the coupling constant u in the vicin- 
ity of d uc — 4 + m/2. In the same way one obtains 
the large- N results at dominant order [32] . More impor- 
tantly for our purpose we recover the low-temperature T 
results obtained in the vicinity of the lower critical di- 
mension di c (m) = 2 + m/2 by Sak and Grest [22] using 
a large- k expansion since K ~ 1/T . We get the flow of k: 
d f K = -2eK + C(N - 2) with C = T[m/A}/2 d TT d / 2 T[m/2] 
and e = d — di c (m) with the running anomalous dimen- 
sions : 77 1 1 = C/k and ?7_l = CzjU. At the fixed point 
one has k* = C{N - 2)/2e and 77^ = 2e/(N - 2) and 
77^ = ??j|/2 which coincide exactly with the expressions 
of Sak and Grest. While fully expected this result is par- 
ticularly valuable for anisotropic systems for which the 
lower critical dimension di c (m = 1) = 2.5 is especially 
close to the physical dimension d = 3. 

We now specialize to the N = 3 case, and thus, m = 1. 
One finds a nontrivial fixed point with two directions of 
instability - corresponding to LCB - in any dimension 
between di c = 2.5 and d uc = 4.5. Fig.l displays the 
curves 77y and 77^ as functions of d, that call for sev- 
eral remarks. First they show the ability of the NPRG 
approach to interpolate smoothly between di c and d uc . 
Second they confirm the salient fact that LCB in d = 3 
is characterized by a negative value of r]\\ a result which 
was often claimed in absence of any direct investigation. 
This result also points out the limits of the pcrturbative, 
large-./V or low-temperature approaches that lead to a 
positive value of 77|| in d = 3. We now focus on the d = 3 
case. For each fixed point the critical exponents isu, u±, 
77|| and 77^ are computed and optimized [28, 33]. To do 
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FIG. 1: The anomalous dimensions m and rj± as functions of 
the dimension d using a field truncation up to 4> s . 

this one considers one (or more) family of cut-off func- 
tions indexed by a real parameter A: (qii)- Typi- 
cally one has considered a cut-off function of the form 
^k ± ( c l\\) = AZ||/(exp(qjj/fcj_) — 1). For each family one 
varies A in order to find stationary values of the critical 
quantities. Stationarity is a condition that must neces- 
sarily be fulfilled by any putative physical quantity to 
ensure its quasi-independence with respect to both the 
cut-off function and the truncation used [28]. However 
an explicit study of the convergence is necessary to get 
trustable results. This has been realized by adding suc- 
cessively powers of the field up to order <fr 12 . Doing this 
we have been able, for all critical exponents, at almost 
any order of the field expansion [34] to find stationary 
values. This is illustrated, for instance in Figs. 2 and 3 
which represent the critical exponents u± and rjx in the 
vicinity of their stationary values for different truncations 
of the action. 
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FIG. 2: The exponent u± as function of A for truncations from 
cf) 4 - lower curve - to (j> 12 - upper curve. Stationary points 
are indicated by black diamonds. 



FIG. 3: The anomalous dimension r]± as function of A for 
truncations from 4 - upper curve - to (j> 12 - lower curve. 
Stationary points are indicated by black diamonds. 

3 also indicate a rapid convergence of the physical quan- 
tities when powers of the field <f> are added. Between 
order <p w and order <p 12 only the third digit of u± and 
rj± are modified (Note however that the parallel expo- 
nents fluctuate a little bit more). This fact indicates 
that our results arc well converged. They are summarized 
in Table (I). An important outcome of our approach is 
that the critical exponents strongly differ from those ob- 
tained within the weak-coupling perturbative approach 
[17]. This discrepancy is not surprising as the perturba- 
tive results have been obtained only at low order. 
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0(0 [17] 




0.78(1) 


0.372 


v± 


1.655(5) 


0.757 




-0.18(2) 


-0.020 


T)± 


0.075(1) 


0.042 



TABLE I: Critical exponents v and r/ for TV = 3 and m = 1. 

To summarize we have shown that the NPRG ap- 
proach provides critical exponents within a systemat- 
ically improvable process that avoids the drastic diffi- 
culties encountered using perturbative approaches. Our 
work could stimulate, in particular, extensive numerical 
works as well as investigation of high-quality magnetic 
compound to confirm the adequacy of our quantitative 
predictions. 

K. E. and D. M. are grateful to B. Delamotte for useful 
discussions. 



Note that the critical exponents vary very smoothly 
with the parameter A around the stationary points which 
indicates a very weak dependence of the results with re- 
spect to the cut-off function used. This fact has been 
confirmed by using other families of cut-off functions that 
lead to the same results. More importantly Figs. 2 and 
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